Looking at
OUTzu_methratioBSMAP_v9_M_A.txt
via
python methratio.py -d /Volumes/NGS\ Drive/Oyster\ Genome/oyster.v9_M.fa -z -u -o /Volumes/AquaculX/armina/OUTzu_methratioBSMAP_v9_M_A.txt -s /Volumes/Bay3/Software/samtools /Volumes/AquaculX/armina/BSMAP_output_trimmed_v9_M.sam
Stats on column 8 (CT count)
| #sum |
mean |
stdev |
0% |
25% |
50% |
75% |
100% |
| 1.46618e+08 |
8.26746 |
23.8425 |
1 |
2 |
3 |
7 |
5929 |
|
original BS seq data
/NGS Drive/NGS Raw Data/Oyster_BSseq/filtered_Unlabeled_NoIndex_L003_R1.fastq
Will Run BSMAP on MBD-BS full genome
hummingbird
./bsmap -a /Volumes/web/cnidarian/filtered_SB_MBD_BSseq_NoIndex_L003_R1.fastq.gz -d /Volumes/web/whale/ce_bs/oyster.v9.fa -o /Volumes/web/cnidarian/BSMAP_MBD_mg_output_v9_A.sam -p 6
Total number of aligned reads: 120734940 (86%)
Done.
Finished at Thu Dec 20 18:10:20 2012